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SYMBOLS 


a vector containing all of the elements of A 

A plant matrix with elements 

b Input distribution vector with elements b. 

C coefficients matrix with elements c^ 

c final-value vector with elements c^ 

C* normalized longitudinal handling-quality criterion 

C*(t) longitudinal handling-quality criterion 

D( ) denominator polynomial 

D* normalized lateral handling-quality criterion 

0* ( t ) lateral handling-quality criterion 

d differential operator d = ^r^- 

G output measurement matrix, matrix transfer function 

augmented output measurement vector matrix 
input/output coupling vector 
augmented input/output coupling vector 
i pilot station to vehicle C.G. distance 

N( ) numerator polynomials vector 

P n normalized roll-rate handling-quality time history 

p(t) roll rate 

q co cross-over dynamic pressure in D* definition 
r(t) yaw rate 

s Laplace transformation variable 

T similarity transformation matrix 

u(t) forcing function 
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V nominal airspeed 

state vector with elements x^(t) 
y output vector with elements y^(t) 

y specified response vector with elements y^(t) 

y augmented output vector 

0 n normalized sideslip handling-quality time history 
0(t) sideslip 

6. aileron deflection 

d 

elevator deflection 
A eigenmatrix 

A specified eigenmatrix 

X eigenvalue vector with elements X^ 

4>( t) roll angle 


INTRODUCTION 


Aircraft control systems have been designed in the past to meet stability, 
apparent damping, and sensitivity criteria. Subsequently, pilot comnents and 
ratings, using the Cooper-Harper rating scale, graded the performance of the 
pllot-alrplane-fllght-control-system combination. The C-H scale emphasizes 

TAKEOFF 

CRUISE 

GROSS MANEUVERING AND AEROBATICS 

FORMATION FLIGHT 

TRACKING 

GROUND-CONTROLLED APPROACH 
LANDING APPROACH [1]* 

Experience has shown that the design criteria are often a subset of the han- 
dling-quality criteria. In order to design flight-control systems which 
achieve low (desirable) C-H ratings, some explicit consideration of handling 
qualities must be Incorporated into the design process. The work begun by the 
Boeing Airplane Company and extended by the McDonnell -Douglas Company has re- 
sulted In quantitative specifications which are Intended to Insure "good" 
aircraft handling qualities. These specifications take the form of envelopes 
within which selected time-histories must fall, e^.g,. , normalized sideslip, 
normalized roll rate, and two blended quantities, C* and D*, containing 
the acceleration cues at the pilot station. 

An evolutionary process has provided the basic design of most aircraft 

control systems. It is often the case that the selection of numerical values 

_ 

The numbers in square brackets refer to bibliography entries. 
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for the many parameters In a control system Is more difficult than the 
establishment of the control system structure and modes of operation. To 
overcome this "tuning" problem, a model -matching optimization technique can 
be Implemented on a large digital computer. Such a program adjusts specified 
parameters of a simulated airplane-flight-control -system combination so as to 
minimize some measure of the dynamic differences between the closed-loop air- 
plane and a low-order model of a desirable prototype airplane. As a partic- 
ipant In the ASEE-NASA Summer Faculty Fellowship Program In the sumners of 
1974 and 1975, the principal Investigator, working in the Vehicle Dynamics 
and Control Division of the NASA Flight Research Center, Edwards, California, 
developed a method of translating handling-quality time-history specifica- 
tions Into prototype aircraft models. These models are suitable for subse- 
quent use with an FRC model -matching program: CONOPT. 

The connection between specified time histories which fall within the 
established envelopes and numerical values for adjustable parameters onboard 
the aircraft Is, in principle , established. The method is explained in the 
following sections and the lateral -motion case demonstrated. 

PROBLEM STATEMENT 

The NASA Flight Research Center computer program CONOPT, which resulted 
from the addition of the MIT Model Performance Index Design Program OPT to the 
in-house control system analysis program CONTROL, is a model-matching package. 
The user specifies the model in transfer-function form. The Model PI tech- 
nique [2] requires that the model transfer function, if it has any zeros, 
should have the same number of excess poles over zeros as the closed-loop 
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airplane control system. If it has no zeros then the number of poles of the 
model transfer function should be equal to or less than the number of excess 
poles over zeros of the closed-loop airplane. Another condition on the model 
transfer function Is that It must have reasonable eigenvalues. Thus the model 
Is a linear, constant-coefficient ordinary differential equation in operato*" 
notation. 


D(d)x(t) = N(d)u{t) , (1) 

with loose bounds on the roots of, 

DU.,) = 0 . (2) 

and constraints on the coefficients of N(d). In order to obtain the entire 
matrix transfer function in one operation the model will be represented in 
state variable form. 


dx = Ax + bu . (3) 

The problem Is to determine an A and b combination which satisfies the loose 
bounds and constraints above and which describes a low-order prototype air- 
plane model with "good" handling qualities. FRC program CONTROL can be used 
to obtain transfer functions from A and b^ for use by CONOPT. 

The handling-quality time histories, y,(t), are linear combinations of 
state variables. 


y(t) - Gx(t) . (4) 

If all the y,{t) fall within their handling-quality time-history envelopes the 
model is appropriate for use with CONOPT. 


f T 


The rigid-body equations of motion for symmetrical airplanes slightly 
disturbed from straight and level flight can be separated Into two sets of 
four simultaneous first-order linear constant-coefficient ordinary differen- 
tial equations [3]. One set describes the longitudinal motion, or motion In 
the plane of symmetry of a normally configured airplane. The other set de- 
scribes the lateral motion. In each case, 

dx = Ax + bu , (5) 

where A Is 4 x 4 and b Is 4 x 1 If there Is only one Input. 

Only one longitudinal handling-quality criterion has been established: 

£ If 

C n< Three lateral handling-quality criteria have been proposed: 0^, normal- 
ized roll rate, and normalized sideslip. An envelope for the first derivative 
of each handling-quality time history has also been proposed [4]. These eight 

envelopes are plotted In Figures 1 through 4. These handling-quality tlme- 

* 

history envelopes imply standard Inputs. The C response results from a step 
change In elevator position, 6 » and the lateral handling-quality time re- 
sponses result from a step change In aileron position, <5 a - 

METHODS 

An obvious procedure for obtaining models with handling-quality time 
histories which fall within specified envelopes Is to: 

1. Guess a model In terms of A and b 

2. Calculate x(t) 

3. Form y(t) 

4. Compare y(t) and y(t) to the appropriate envelopes 
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5. If no envelopes are violated, stop 

6. Otherwise alter A and/or b according to some 
strategy and return to Step 2. 

This "brute force" technique would consume large amounts of computer time and 
may not succeed. It would be necessary to define a cost functional which 
measures all excursions outside of the envelopes and which is suitable for use 
in a numerical optimization scheme. Next a connection between the elements 
a^ and b^ and the cost functional would have to be established. Finally, 
there is no assurance that the attainable minima of this functional would be 
zero. The addition of hard and soft constraints arising from the loose bounds 
on the eigenvalues and specification of the pole-zero excess of the transfer 
functions complicates the numerical optimization problem. In the fourth-order 
lateral -motion case there would be: twelve variables to be manipulated, four 
soft constraint equations to be approximately satisfied, n hard constraint 
equations where n is the sum of the specified pole-zero excesses, and six en- 
velopes to be matched. 

One simplification which circumvents the need for a cost functional and 
optimization strategy is the use of direct search. If any solution which sat- 
isfies the constraints and the envelopes Is as useful as any other, systematic 
or random direct search is easier to Implement and no more wasteful of computer 
time than optimization. Random direct search has been used to obtain second- 
order longitudinal airplane transfer functions. In this simpler situation 
the a^ : i,j = l,2 are randomly selected. If the resulting eigenvalues fall 
within the region shown in Figure 5 [4], the b^ : i = 1 ,2 are randomly selected 
and c* and C* calculated. The ranges of the random choices are very 
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Influential. Sample results are: 

1. 558 random sets of a^ yielded 22 acceptable pairs of real 
eigenvalues and six acceptable pairs of complex eigenvalues. 

2. For A having acceptable real eigenvalues 118 random pairs of. 
b,| were required to find one pair which yielded time-responses 
which match the C and C envelopes. 

3. For A having acceptable complex eigenvalues 320 random pairs 
of bj were required to find one pair which satisfied the 
envelopes. 

4. If one envelope was matched (and the other violated) it was 
as likely to be the C envelope as the C envelope. 

The application of direct search to the fourth-order problem is currently 

being Investigated. This is potentially useful for the development of models 

* * * 

which satisfy the longitudinal criteria, C and C , but less so for the lateral 
ease. In order to develop lateral motion models and obviate the use of iter- 
ative numerical techniques the problem must be reposed. 

If one begins with specific time responses which fall within the envelopes 
and seeks a model whose output responses closely resemble the specified 
handling-quality time histories the problem no longer requires iteration. If 
the eigenvalues are specified or obtained from the input responses and then 
frozen and if all elements of A and b are unspecified rather than just those 
elements which are normally nonzero for most airplanes, the problem is well- 
posed. This procedure has three disadvantages: 

1. The partial degrees of freedom represented by loose bounds on 
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on the eigenvalues are lost. 

2. The envelopes are not directly used and may be violated by 
the results. 

3. The resulting A is not constrained to have the zero elements 
and kinematic elements usually found In such models of 
airplanes. 


Experience with lateral -motion examples has shown that A may well have ob- 
viously unrealistic elements and still produce realistic transfer functions. 
This Is sufficient for present purposes. 

The reposed problem Is the following: given y(t), h(a^j) and G(a^ j ) , 
find A and b such that if 

a* m 


and If 

then 

with 


dx = Ax + bfi a , 
y = Gx + h<$ a , 


A 

y s y In some best sense, 

A 

A - A 


( 6 ) 

(7) 


( 8 ) 


Note that the elements of G and h are known functions of the unknown elements 

of A and that G Is normally not square, there being fewer elements in y than 

in x. 

«* 

The first step is to obtain an analytical representation of the input 
time histories which one specifies in the form of a table of discrete values 
uniformly spaced in time. If the model is to be a fourth-order constant- 
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coefficient linear ordinary differential equation. Its solution must have 
the form: 

y^t) = c i ie Xlt + c 12 e Xzt + c i3 c X3t + c i4 e X4t + c 1 . (9) 

These coefficients and eigenvalues, and Aj, are obtained by least-squared- 
error fitting the above form to the specified discrete values of y^. This is 
a two-step process. First the Aj are calculated. The calculation is based 
upon first and second differences between adjacent y^ values and Is strongly 
Influenced by their precision. If the values are represented by three signif- 
icant figures the resulting eigenvalues are drastically incorrect even when 
the data are contrived and there should be an exact solution. As the number 
of significant figures is increased the calculated eigenvalues more closely 
resemble the correct values. However, if data obtained from sketched time 
responses or Imprecise tables are to be accommodated, it is impractical to 
calculate eigenvalues from such data. In such cases, the eigenvalues which 
the model is to have must be specified. In either case, the c^ are straight- 
forwardly calculated and the resulting fits usually quite good. Thus one can 
represent the Input time histories, y(t)» as 

y * y = Ce Xt + c (10) 

where the elements of c are the c^ . 

An example set of discrete input data points and their fitted represent- 
ations Is shown In Figures 6 through 8. Reasonable responses were sketched 
onto graphs of the lateral handling-quality envelopes and then converted to 
tabular form by estimating the values at half-second intervals. The eigen- 
values used In the curve-fitting process were specified. The values used were: 
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TIME SECONDS 

FIG. 6 FITTED ROLLRRTE RESPONSE 
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\i - -2.4 

X 2 = -.003 

X 3 * -.25 - j2.0 

X„ = -.25 + j2.0 

A 

Once the specified discrete values of the are approximately represented 

y = Ce-* + c 

This equation can be Laplace transformed Into 

!<*> -frlr • on 

Assuming a unit-step Input, the plant and output eq nations can also be Laplace 
transformed if to 


sX(s) = A X( s ) + b/s 
Y ( s ) = GX(s) + h/s 
which are manipulated to obtain 

X ( s ) = (si - A)“ l b/s = $(s)b/s 

Y(s) = G<£(s)b/s + h/s . (12) 

If equation 12 is rearranged to match the form of equation 11, the right- 
hand sides of equations 11 and 12 can be equated, coefficient by coefficient, 
to produce n(m+-1) simultaneous algebraic equations where n is the order of 
the plant equation, equation 6, and m is the number of elements in the output 
vector, equation 7. For low-order models this is a satisfactory technique. 
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The algebraic equations are of degree <n. In the fourth-order lateral - 
motion example there are sixteen nonlinear equations containing, at worst, 
quadruple-cross-product terms. A limited amount of experience with these 
equations has shown them to be numerically sensitive. An effort Is being 
made to develop a computer program which will solve the fourth-order example 
equations with various constraints on the elements and eigenvalues of A. 

While this may prove to be the most satisfactory approach to the problem, the 
numerical uncertainties are sufficient to motivate the development of an al- 
ternative technique which Is computationally simpler. 

A less numerlcally-sensltlve approach can be developed by combining equa 
tlons 7 and 10. 


y = Gx + h6 a 
x = G' l y - G _1 h6 a 
x = G^Ce^ + G -1 c - G'‘ 1 h6 fl 
dx = G^CAe-* - G _i hd$ a 

If these expressions are substituted into equation 6 one obtains 

G“ 1 CAe- t - G _1 hd<S = AG^Ce- 1 + AG _1 c - AG^hfi^ + b6 3 
or 

(G -1 CA - A6"*C).e^ t = AG _1 c + G -1 hd6 a + (b - AG’ 1 h)6 a . 

For this expression to be true 

G"*CA = AG _1 C (13a) 


20 


(AG” l h - b)6 = AG“‘c 

•• • - a •• ■» 


(13b) 


If I » G _1 C, then 


and 


A * TAT -1 


(14) 


b - AG" l (h - c) . (15) 

It Is assumed that 6 a (t) = 1 and G^hdS. Is neglected at this point but must 
be introduced as an Initial condition vector when computing model time re- 
sponses. A model obtained from equations 14 and 15 will have the specified 
eigenvalues and match the analytical representations of the Input data exactly. 

Unfortunately the above results assume the existence of G -1 whereas, in 
general, G is not square. One can substitute the right pseudoinverse and ob- 
tain similar results except that T will be singular and the model eigenvalues 
are no longer equal to the specified eigenvalues since A and A are no longer 
similar. If G Is square and nonsingular, equations 14 and 15 hold. There are 
two cases of practical Interest in which G is nonsingular. First, one can 
soeclfy n 3 m, that is, let the number of time histories establish the order 
of the model. Secondly, it may be possible to adjoin created pseudo-time- 
histories to y until m 3 n. One must have the means to create these pseudo- 
time-histories. They must be linear combinations of the model state varia- 
bles and the resulting distribution matrix must be nonsingular. 

In the lateral -direction, small -motion example, the state variables are 
roll rate, yaw rate, sideslip angle and roll angle. The handling-quality 
Indicators are normalized roll rate, normalized sideslip angle and normal i zed D*. 
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The normalized roll -rate Input data can be Integrated to produce normalized 
roll -angle pseudodata. For the example, this was accomplished by exactly 
fitting a tenth-order polynomial to the eleven discrete roll -rate data points 
shown In Figure 6. This polynomial was Integrated and evaluated at the same 
values of time to generate the discrete pseudodata points shown in Figure 9. 

The continuous curve In Figure 9 results from fitting an equation of the form 
of equation 9 to the discrete pseudodata points just as is done for the normal 
data. This allows the expansion of y to y and G to G and assures the pon- 
singularlty of G. This, In turn, allows calculation of A and b using equa- 
tions 14 and 15. The same straightforward opportunity to augment y until G 
Is square and nonsingular does not exist in the longitudinal small-motion case, 
although a designer could Impose specifications on the motion, in addition to 
the Boeing C* criterion, and thereby create a nonsingular G. 


COMPUTATIONAL PROCEDURES 

D* is a weighted combination of aircraft sideslip which is considered 
the principal low-speed handling-quality parameter and lateral acceleration 
at the pilot station which is the principal motion cue parameter at high 
speeds [4]. In terms of a stability-axis system representation of small 
perturbations about straight and level flight for a normally configured air- 
craft with the pilot station on the longitudinal principal axis the expression 
for D* is 

D*(t) = V(B + r) + ir + c 3 q G0 B 06) 

's 

where Is a dimensional constant defined in Table 1. In terms of the 
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Note that this assumes that the plant matrix Is 


A a 


a ll a 12 
a 21 


With no constraints on the elements a^. Thus A Is specifically not con- 
strained to be of the form 


similarly 


A = 


L p L r ® 
N p N r N & 0 

a 0 -1 Yfi Y 


1 o o 


(19) 


b = 


and not 


Thus the output distribution matrix 6 for the lateral motion case is of the 
form 


G = 


1/P 


ss 


V D ss 


0 

0 

V D ss 


0 

1/Bss 

V D ss 


0 

0 

Da/D, 


ss 


( 20 ) 


where 


y 1 (t) = p(t)/p ss = P n (t) 


y 2 (t) 


3(t)/e ss - 0 n (t) 


origins 
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y 3 (t) = D*(t)/D ss = D*(t) 


( 21 ) 


and 


h = 


The augmented output distribution matrix is 


0 

0 

V D « 


( 22 ) 


ss 


1/p 

0 

0 


0 / D 
p' ss 


0 0 0 

0 1/3 SS 0 

0 0 Vp ss 

V D ss V D ss V D ss 


(23) 


which is nonsingular if D r f 0. The augmented output vector, y(t), is 


y-,(t) = p n (t) 
y 2 (t) = e n (t) 

(24) 

y 3 (t) = | p n (t)dt = d> n ( t) 
y 4 (t) - D*(t) 

o 

Since n elements of A and the n elements of b are to be established equation 

2 2 

14 Is equivalent to n scalar equations with n unknowns and equation 15 is 
equivalent to n scalar equations with n unknowns. This is not significantly 
different from the Laplace transform method in which one obtains n(m+l) alge- 
braic equations without having had to create psuedodata. However, the degree 
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of the simultaneous equations to be solved In the Laplace transform method 
Is n whereas the equations arising from the pseudodata method contain cross- 
product terms at worst. This dramatically Increases the llkllhood of obtain- 
ing solutions by iterative numerical means. For the fourth-order lateral 
motion example, there are sixteen equations to be solved for the a^. A 
Newton-Euler procedure [5] was employed in which an initial estimate of the 
solution vector a Is Iteratively improved. If equation 14 is rewritten 


f(a) = 0 


(25) 


and P(a) is the Jacobian matrix associated with the equation 25 then 


2m ■ 2 k - rVidi.) 


(26) 


This simple procedure has proved to be so successful that A = I can be used as 

-5 

the starting point and every element of (a^ - a^) Is reduced to 10 in three 
iterations typically. The elements of b are obtained from equation 15 without 
iteration. 


ft’jfore A and b can be calculated, equations of the form of equation 10 
must be fitted to the input data. The resulting coefficient arrays C and c 
appear in equations 1*> and 15. Note that T * E C in the pseudodata method. 

A set of eigenvalues can be obtained from each discretized input time 
history. If 


e z = y^zAt) - y^(zAt) 


and 


= e X J At 


original pagp tc 
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then 

e z ° C 11 u 1 + C 12 u 2 + C i3 lJ 3 + C 14 y 4 + C i ' (27) 

where 

z = 1,2,... (number of discrete values y^(zAt)). 

The c^j and can be eliminated by linearly combining the e z « One obtains a 
linear set of simultaneous algebraic equations 

D’y * d‘ (28) 

where the elements of D’ are first differences of the table of y^ (zAt) values 
and the elements of d' are second differences. Unfortunately, if the y^(zAt) 
values are obtained from time-response sketches or are imprecise for any reason 
the p. calculated from equation 28 are uselsss. Frequently the Uj obtained in 

j J 

such cases have negative values from which no X. can be calculated. 

J 

In lieu of eiganvalues obtained from the input data, it has been neces- 
sary , in all practical calculations, to use specified eigenvalues. In such 

cases the calculation of C is based on minimizing 

i i 

Z" i 

where the lateral handling-quality time histories typically span five seconds 
and are discretized by taking values every half second for a total of eleven 
values per time history. The ir :1al value of y(t) is forced to be zero by 

C 1 ' j^j C i j * (t29> 

As a check on the entire computational process once an A and b have been 
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calculated, G and h can be calculated, x(t)obta1ned by Integration and y(t) 
calculated. This y{t), obtained by Integration, Is plotted with the fitted 
y(t). They should be identical. The Integration method is described in 
Takahashl, et. al_. , page 103 [6]. The solution matrix Is represented by a 
series expansion containing p terms. The number of terms is specified by a 
recipe attributed to Paynter 


~r (nq) P e ng 2 0.001 {30) 

where n is the size of the A matrix and q is the largest element of AAt. 

The example fitted and integrated handling-quality time histories are 
plotted in Figures 10, 11, and 12. 
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ZERO SUPPRESSION 

The Laplace -transformed analytical representation of the Input data, 
equation II, has one fewer zero than poles in each element of the matrix 
transfer function G( s) where 


Y(s) a G(s)/s 


or 




This expression can be rewritten, omitting the 1 subscript 


( C 1 X 1 +C 2 X 2 + C 3 X 3 + C 4V S “ ( C 3 X 3 X 4 + C 4 X 3 X 4 + cX l X 3 


^ X^3 + c i a i c 4 X 1 ^ ^*2 X 2 X 3 ^ ^*2 X 2 X 4 ^ ^4 X 2 X 4 ^ 1 X 2 X 1 


* c 2 * C^A-j * C^X^XqA 


1 A l A 3 4 3 1 3 4 4 1 3 4 u 4 1 3 4 C 2 A 2 3 4 


+ c 3 A 2 X 3 X 4 + c 2 X 2 X 3 X 4 + C 1 X 1 X 2 X 3 + C 2 X 1 X 2 X 3 + C 3 X 1 X 2 X 3 


+ c 


1 X 1 X 2 X 4 C 2 X ] X 2 X 4 c 4 X i X 2 X 4 ) ® ^ ^7 X 2 X 3 X 4 


(31) 


To reduce the number of zeros in the transfer function, the following con- 
straint equations must be incorporated Into the least-square minimization, 
equation 29. To suppress one zero 


c i x l + C 2 X 2 + C 3 X 3 + C 4 X 4 0 * 


(32) 


To suppress two zeros 


c 7 A i + C 2 X 2 + C 3 X 3 + C 4 X 4 0 


and 

'34 
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C 3 X 3 X 4 + C 4 X 3 X 4 + C 1 X 1 X 3 + C 3 X 1 X 3 + C 1 X 1 X 4 + C 4 X 1 X 4 

+ C 2 X 2 X 3 + C 3 X 2 X 3 + C 2 X 2^4 + C 4 X 2 X 4 + C I X 1 X 2 + C 2 X 1 X 2 * 0 ( 33 ) 

If three zeros are suppressed, only one parameter remains to be established 
by the minimization and the resulting representation of the input time history 
Is inadequate. 

For low-order models, zero suppression has two detrimental effects on the 
models obtained. The first is the worsening of the agreement between the dis- 
cretized input time histories and their analytical representations. For ex- 
ample, the unconstrained least-squa red-error method yields a satisfactory fit 
of the roll-rate Input, as shown in Figure 6. The continuous fitted curve has 
approximately the same relationship to the envelope as does the discretized 
Input data. If the LSE method is constrained by equation 32 to have one less 
zero, one obtains the fit shown in Figure 13. If one further constrains the 
LSE minimization by equation 33, another zero will be eliminated and, for the 
example data, the fit Is degraded to that shown In Figure 14. The same se- 
quence is portrayed in Figures 7, 15, and 16 for the sideslip input. The 
Figure 7 results are unconstrained, one zero is suppressed in Figure 15 and 
two zeros are suppressed in Figure 16. Important features of the discrete 
Inputs can be preserved by weighting the appropriate errors, e z , in the LSE 
method. This may make the constrained fits more useful. Even the uncon- 
strained cases can be altered. The risetime of the fitted roll rate can be 
reduced, for example. However, the overall fit cannot be improved by weighting. 

The second detrimental effect of zero suppression is the increased ten- 
dency of the A matrix to contain unrealistic values. The A and b arrays re- 
sulting from unconstrained fits of the input time histories are 
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-1.36 

-4.17 

-20.22 

-.02 

-1.12 

1.89 

12.43 

-.04 

0.22 

-.80 

-3.43 

0.02 

0.20 

1.04 

5.61 

0.00 


The Newton-Euler procedure far obtaining the converged in three iterations 
to a maximum variation of any element from one iteration to the next of 
.00001 or less. This model was obtained from the unconstrained fits shown in 
Figures 6, 7, and 8. 

If one zero in the roll -rate transfer function and one zero in the side- 
slip transfer function are suppressed, the A and b arrays are altered to 


r 

- 3.49 

4.73 

97.14 

0.52 

0.96 

- 3.42 

- 22.48 

-.07 

-.19 

0.30 

4.00 

0.02 

-.32 

2.44 

15.42 

0.01 


The Newton-Euler procedure again required three iterations. This model was 
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obtained from the fits shown In Figures 13, 15, and 8. 

If two zeros In the roll -rate transfer function and two zeros In the 
sideslip transfer function are suppressed, the A and b arrays are 


-9.34 

-1.66 

403.88 

2.75 

3.51 

-2.88 

138.62 

-.86 

-.21 

-.04 

8.73 

0.06 

-2.10 

1.45 

96.10 

0.59 


and 

b = 

Three Iterations were sufficient and the model was obtained from the data 
shown in Figures 14, 16, and 8. 

The number of terms in the expansion of the state transition matrix speci- 
fied by Paynter's recipe, equation 30, is a function of the largest absolute 
value contained in A and this number is frequently larger than the practical 
limit of about 30 when zeros are suppressed. 

The additional information required to calculate A and b from the discrete 
time histories shown in Figures 7-16 is given in Table 2. The specified 
eigenvalues were 

A 1 = ~ 2A ORIGINAL PAGE IS 

OF POOR QUALITY 

X 2 = -.003 
X 3>4 = -.25 + j 2 . 0 


0.00 

1.45 

0.00 

0.87 
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AIRPLAKE MODEL WITH SPECIFIED TIME HISTORIES 


I 



.42 


i 


l 


The eigenvalues and vehicle and flight parameters are approximately those of 
a Lockheed Jetstar, a four-engined utility transport [7], at 20,000 feet al- 
titude and Mach = 0.6. It should be noted that the discrete input time his- 
tories were obtained fromp n , and D* responses sketched by an FRC engineer 
[8] and do not refer to a particular airplane or flight condition. Finally, 
FRC program CONTROL was used to calculate transfer function coefficients from 
the above A and b arrays In combination with the appropriate distribution 
matrices, calculated from equation 22, to verify the suppression of the speci- 
fied zeros. 


VERIFICATION 

Flight test data obtained from the FRC Lockheed Jetstar [9] provides a 

test case for the above model -generation procedure. For small, lateral varia- 

★ 

tions about straight and level flight, the Jetstar can be represented by 


-2.353 

0.735 

-11.050 

0.000 

1 

O 

tn 

0.358 

3.836 

0.000 

0.026 

-.999 

-.205 

0.053 

1.000 

0.054 

0.000 

0.000 


5.650 
0.031 
-.001 
0.000 

Body axis nondimensional stability derivative parameters used in place of 
stability axis values for verification purposes only. 
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1. 

0. 

0. 

0. 

0. 

0. 

1. 

0. 

0. 

0. 

0. 

1. 

14.75 

-7.044 

-146.0 

32.14 


L. 


— 

0 . 

o. 

- = 0 . 

0. 

The eigenvalues were 

X 1 = -.24045 
\ 2 = -.00310 

X 3 4 = -.25428 + j2. 06475 


CONTROL was used to calculate the handling-quality time histories for the 
period 0.0 to 5.0 seconds. The discrete input data extracted from the CONTROL 
output is given in Table 3. The flight and vehicle parameters are given in 
Table 2. The normalization factors were arbitrarily selected. 


The responses of the resulting model are shown in Figures 17, 18, and 20 


and compared with the discrete input values of Table 3 after normalization. 


The A, b, G and h arrays which define the model are 


-2.359 

0.771 

-10.939 

-.003 

-.053 

-.359 

3.849 

0.000 

0.025 

-1.005 

-.201 

0.053 

0.906 

-.021 

.278 

0.003 


A 
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I'I'l 




5.818 


b = 


-.053 

-.215 


0.050 




L_ 

_i 



1.0 

0. 

0. 

0. 


0. 

0. 

1. 

0. 

G ** 

0. 

0. 

0. 

1. 


14.35 

-11.13 

-143.5 

32.61 



These results are close to the original arrays. The differences are attrib 
uted to the truncation of the CONTROL output values to two or three signifi 
cant figures to approximate imprecise or sketched input data. 
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DISCUSSION 


The problem of obtaining an A and b from specified output time histories 
Is one of ,*>nl inear, noisefree identification. Five techniques for solving 
this problem have been suggested. The minimization of a cost functional which 
measures the differences between a trial solution and the handling-quality 
time-history envelopes would consume a large amount of computer time and there 
is no assurance that such a cost functional would be sufficiently well behaved 
to have a useful solution. Similarly, one could consume a large amount of 
computer time seeking solutions by random direct search. A graduate student 
is currently working on a variation of direct search in which the sensitivity 
of the time-history errors to changes in the A matrix elements is calculated. 
Then a set of incremental changes to the A elements can be obtained. 

The Laplace transformation method is also being pursued by a graduate 
student. In this form the problem is quite easily formulated but is i'il- 
posed. It remains to be seen whether this method can advantageously incor- 
porate the loose bounds on the eigenvalues or not. It also has the disad- 
vantage of being a two-approximate-step process. One obtains y ( t ) from y(t) 
and then y(t) from y(t) except that the latter two will not coincide as they 
do in the pseudodata method. The Laplace transformation method and the sensi- 
tivity matrix method do have the advantage that they can be made to yield A 
matrices of the form of equation 19. 

The remaining approach, the pseudodata method has two advantages. It 
contains only one approximation step and it is numerically efficient. The 
disadvantages are that it is somewhat less general and yields unconstrained A 
matrices. Experience has shown that the transfer functions resulting from 
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the unconstrained A matrices resemble those produced by A matrices of the 
form of equation 19. The pseudodata method results in a well -posed set of 
bilinear algebraic equations which yield an A matrix having the specified 
eigenvalues. 

The pseudodata method is the only method which readily achieves zero 
suppression. In the other methods zero suppression contributes to their i 1 1 - 
posedness making useful solutions even more numerically difficult to obtain. 
The utility of any zero-suppressed solution is called into doubt by the 
detrimental effect suppression has on the LSE fit of the discretized input 
data (see Figures 13, 14, 15, and 16). Relieved of the need to suppress 
transfer-function zeros, one might prefer one of the other methods. 
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APPENDIX A 

Use of Program AANDB 


INTENDED OPERATION 


t 


t 


1. Put in normalized discrete data 

2. Obtain normalized responses 

3. Put in normal ization factors 

4. Obtain the normalizing distribution matrix 

5. Obtain the non-normal ized airplane equations 
in state-space form. 

6. Use CONTROL to obtain transfer functions 

Steps 1 through 5 are performed by an example Fortran program AANDB for 
the fourth-order small -lateral -motion case. The structure of the program is 
shown in Figure A-l . The subroutines perform the following tasks: 

MUGEN: Calculates eigenvalues from input time-history data. 

CASE! * (with entry points CASE2 and CASE3) : Calculates the c^ 

required to fit c^e^ 1 ^ + Cpe^ 2 ^ + c^e^ + c 4 e ^ t + c 5 
to the input time-history data. 

LISTER: Prints information from COMMON on demand. 

FITTING: Fits a polynomial to roll-rate data and integrates it to 

produce roll-angle pseudodata. 

MODEL: Calculates the plant matrix, A, input distribution vector, 

b, the output distribution matrix, G, and the input/output 
coupling vector, h_. 

RESPONS: Integrates the plant equations to produce comparison time 

histories. 

CASE1 does not suppress any transfer function zeros, CASE2 suppresses one zero 
and CASE 3 suppresses two zeros. 




READ- 


READ -- 




MUGfN 


FITTING | 

— 

INVR 





-jGRAPHS PLOT 


LISTER PRINT 





— J 

1 

RESPONS 

LISTER 



1 1 CTCD 



Lio 1 1 K 

— 

MODEL 


INVR 



PRINT 


FIG. A1 
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GRAPHS: Plots envelopes and time histories and their first derivatives. 


1NVR: Inverts matrices. 


The first fourteen input data cards are read by MAIN. The final two 
input data cards are read by subroutine MODEL. The input data cards contain 
the following information: 


1st data card : 
IF NN(1) = 1: 

NN(2) = 1: 

NN(3) = 1: 
NN(4) = 1: 

NN (5) = 1: 
NN(6) = 1: 
NN(7) = 1: 
NN(8) = 1: 
NN(9) = 1: 

NN(10) = 1: 

NN (11) = 1: 

NN (12) = 1: 
NN( 1 3) = 1: 
NN (14) = 1: 

NN (15) = 1: 


NN( 20) ; 2012 

LISTER will print: data, NN, MCASE, NPLOT, NHIST , 

input time histories 

LISTER will print eigenvalues obtained from input 
time histories 

LISTER will print eigenvalues specified by user 

LISTER will print roll angle data obtained from 
roll rate 

LISTER will print complex coefficient matrix C 

LISTER will print summary of time-history curve fitting results 

not used 

not used 

LISTER will print envelopes for PN, BETAN, DSTAR 
& derivatives 

LISTER will print fitted curves for PN, BETAN, DSTAR 
& derivatives 

LISTER will print number of terms included in series 
for 

LISTER will print difference equation P matrix 

LISTER will print difference equation g vector 

LISTER will print responses obtained by integrating 
model equations 

LISTER will print derivatives of responses obtained 
from P, q 
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NN(16) = 1: 
NN( 1 7) = It 
NN (18) = 1: 
NN (19) = 1 : 
NN(20) = 1: 


LISTER will print "PN PLOTTED" if NPLOT requests it 
LISTER will print "BETAN PLOTTED" if NPLOT requests it 
LISTER will print "DSTAR PLOTTED" if NPLOT requests it 
not used 
not used 


2nd - 12th data cards : TIME, PN, BETAN, DSTAR; 4F10.0 


0.0 

0.0 

0.0 

0.0 

0.5 

_ PN(.5) _ 

BETAN (.5) 

DSTAR( . 5) 

1.0 

PN(1.) 

BETAN (1.) 

DSTAR (1.) 

1.5 

_ PN (1,5) 

BETAN(1 .5) 

DSTARfl .5) 

2.0 

PN(2. ) 

BETAN (2.) 

DSTAR (2.) 

2.5 

PN (2.5) 

BETAN(2.5) 

DSTAR(2.5) 

3.0 

PN{3. ) 

BETAN (3.) 

DSTAR{3. ) 

3.5 

PN(3.5) 

BETAN(3. 5) 

DSTAR(3.5) 

4.0 

PN(4. ) 

BETAN (4.) 

DSTAR(4. ) 

4.5 

PN (4.5) 

8ETAN(4. 5) 

DSTAR(4. 5) 

5.0 

PN(5.) 

BETAN (5.) 

DSTAR{5, ) 


13th data card : MCASE, NPLOT, NHIST; 314 

If NHIST f 0 subroutine RESPONS will be called to calculate the 
model time histories. If NPLOT f 0 subroutine graphs will be called 
to plot the model time histories according to: 


OKlG^ iVL 

qt? poos 


PAG ft* 



NPLOT 

0 no plots 

1 PN and PNOOT plotted 

2 BETAN and BETANDOT plotted 

3 PN, PNDOT, BETAN, BETANDOT plotted 

4 PSTAR, DSTARDOT plotted 

5 BETAN, BETANDOT, DSTAR, DSTARDOT plotted 

6 PN, PNDOT, DSTAR, DSTARDOT plotted 

7 all plots 

8 no plots 


MCASE must be appropriate for: 


GO T0(1, 2, 3, 4, 5, 6, 7,8,9), MCASE 


1 Pjs) has no zero(s) suppressed 
0 { s ) has no zero(s) suppressed 

2 P(s) has two zero(s) suppressed 

B(s) has one zero(s) suppressed 

3 P(s) has one zero(s) suppressed 

0 ( s ) has two zero(s) suppressed 

4 P ( s ) has one zero(s) suppressed 

B(s) has one zero(s) suppressed 

5 P ( s ) has one zero(?) suppressed 

B(s) has no zero(s) suppressed 

6 P(s) has no zero(s) suppressed 

B(s) has one zero(s) suppressed 

7 P(s) has two zero{s) suppressed 

B(s) has no zero(s) suppressed 

8 P ( s ) has no zero(s) suppressed 

B(s) has two zero(s) suppressed 

9 P(s) has two zero(s) suppressed 

B(s) has two zero(s) suppressed 
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14th data card : specified eigenvalues; 8F10.0 


1-10 First eigenvalues which is real 
21-30 Second eigenvalue which is also real 
41-50 Real part of third eigenvalue 
51-60 Imaginary part of third eigenvalue 
61-70 Real part of fourth eigenvalue 
71-80 Imaginary part of fourth eigenvalue 


15th card: 


velocity, length, c 3 , q cQ , P ss , 




7F10.0 


VELOCITY 

LENGTH 



P 

8 

D 


ss 

ss 

* 

ss 


Nominal airspeed in ft/sec 

Centerline length from CG to pilot in ft 

, * 

Dimensional constant for 0 

2 

Nominal dynamic pressure in lb/ ft 
Roll rate normalization factor 
Sideslip normalization factor 
DSTAR normalization factor 


16th data card : Newton-Euler parameters; 14, F20.0 


ITMAX, 14, Maximum number of iterations of Newton-Euler algorithm 

EPSI, F20.0, when every unknown (the elements of the A matrix) changes 
by an amount smaller than EPSI, the Newton-Euler Algorithm 
stops 

Note: 50,0.00001 seem to work well. 


INTERPRETATION OF THE PRINTED OUTPUT 


If NN - 20*1 all of the following output will be produced: 





LISTER calls from the main program print: 

1. Date, 11s$ table NN, case number, plot request code, time response 
code and the Input time history data. 

2. Real and Imaginary parts of each MU(4) and El (4) obtained from the 

0 5A 

Input time histories. The El are eigenvalues, X^, and each MU Is s e i 
where 0.5 Is the uniform time interval between Input time-history data points. 

3. Real and imaginary parts of the specified eigenvalues and associated 
MU values. These ere the eigenvalues used in the least-square fitting process 
to obtain C, not the quantities derived from the data. The derived values are 
presented only for comparison purposes. 

4. The PHI or roll-rate psudodata generated from the PN or roll rate 
input time history. The time intervals are the same as for the original data. 

5. The TIME RESPONSE COEFFICIENTS matrix, C ’ , is printed. This is a 
4x5 array of complex numbers; 

X = Ce-t + c , where C' = [C*c] , 

which analytically represents the input time histories. This y(t) is printed 
as: FITTED TIME RESPONSES by a call to LISTER from RESPONS. 

6. A summary printout of MU, El, C and c 
LISTER calls from subroutine RESPONS print: 

★ d d 

1. The envelopes, upper and lower boundaries, for P n ® n ’ ^ ’ it p n’ ~5t ^n’ 
and D*. Values are given for every 0.1 seconds. These values are stored 
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in DATA statements in the main program. 


2. The curves which have been fitted to the input time-history data are 

•k t 

printed for p^, $ n and D . The analytic expressions for the curves are differ- 
entiated and tabulated also giving p n , & n and D . 

3. The number of terms taken to calculate the truncated stries used to 
represent e- At is PAYNTERS RECIPE NUMBER. See CONTROL, Takahashi , et. a].. , 
page 103 [6]. 

4. The difference equation parameters P and q in: 

-k+l 3 f -h + 3 u k • -o a 

5. The numerically integrated response time histories y^ where y^ = 

--k + - u k' 

6. The first derivatives of the integrated responses are tabulated at 
0.1 second intervals. They are obtained from 

*k - Ax k + b U|< 
and 

ik = 

LISTER calls from subroutine GRAPHS prints: 

1 If HP l Of ir. such that plots are requested, LISTER will print "PN 


0 

D 6 a /D r 

0 

0 


ORIGINAL PAGE IS 
oi P °QR QUALITY 
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PLOTTED" if p n and p plots have been generated, "BETAN PLOTTED" if B n and 
^ p plots have been generated and "DSTAR PLOTTED" if D* and D plots have 
been generated. The envelopes are automatically added to each plot. 

SAMPLE PLOTTED OUTPUT 

Plots are produced in three groups which can be requested individually or 
in any combination. These groups, the PN group, the BETAN group and the DSTAR 
group, consist of a title line, a lower graph and an upper graph. Each lower 
graph shows the upper envelope boundary, the lower envelope boundary, the 
analytical curve which has been fitted to the input data and the time history 
obtained by integrating the equations of motion. These four traces are repre- 
sented by continuous lines and are plotted versus time on the horizontal scale. 
The fitted and integrated lines should coincide. In addition, the lower graph 
contains eleven discrete symbols representing the input data. 

Each upper graph shows the upper envelope boundary, the lower envelope 
boundary, the first derivative of the fitted analytical curve and the first 
derivative of the integrated time response. The last two should coincide. All 
four curves are represented by continuous lines and are plotted versus time on 
the horizontal scale. A PN group sample is shown in Figure A-2, a BETAN group 
sample is shown in Figure A-3 and a DSTAR group is shown in Figure A-4. The 
example plots are unusual in that the original data was generated by a linear 
simulation program using a fourth-order model. Thus one should expect excel- 
lent agreement between the input data and the fitted and integrated results. 
Data obtained from other sources will not be matched as well. 


-62 



Ill I 









EXPERIENCE WITH THE PROGRAM 


Several bounds on program flexibility exist. The first involves the 
generation of eigenvalues from input-time history data. Subroutine MUGEN does 
calculate a set of A., from each set of eleven points representing one input 
time history. These values are not correct or consistent even when the data 
is carefully contrived for best results. This is due to the dependence of the 
A.j calculation upon differences between adjacent (in time) data values. In 
anticipation of the program being normally used with data obtained graphically, 
the values used in the example were held to three significant figures. This 
so degraded the A^ calculation that the values from the three input time his- 
tories varied from the known values used to generate the data by two orders of 
magnitude in the case of small eigenvalues and occasionally had the opposite 
sign. In general, the complex eigenvalues were more closely identified than 
the real eigenvalues. This was particularly true if the data precision was 
increased. In view of the anticipated character of the input data and the de- 
sirability of specifying the eigenvalues of the resulting model, no further 
development of MUGEN is planned. 

A second difficulty arises when zero suppression is specified. Subroutine 
CASE2 and CASE3 calculate coefficients, C, which yield transfer functions 
having the correct pole-zero excesses. The resulting expression will normally 
be a poor representation of the input time-history data but that is unavoidable. 
Unfortunately, the suppression of zeros also tends to ill -condition the nu- 
merical equations which must be solved to obtain A. In extreme cases this pre- 
vents convergence of the simple Newton-Euler-Raphson algorithm employed in 


MODEL. 
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Finally, there is a general difficulty that may interfere with the in- 
tegration of the model equations of motion in RESPONS. If the A matrix i c . 
ill-conditioned an unworkable number of terms in the series aooroximation for 
e- At may be required. This is caused by the fixed At of 0.1 seconds which is 
built into RESPONS. A limit of 30 terms is imposed. Other limitations will 
undoubtedly come to light as experience with the program accumulates. 


APPENDIX B 


Program AANDB Listing 
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CALL QLINE(TIHE*0STARlL»51,m G*A 7A0 

CALL QLINE(TIME,DSTAR1F,51, 01 G?A 750 

CALL 0LINE(TIP£t0STARlC«51, 01 GR A 769 

CALL ENOPLT GR A 770 
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APPENDIX C 

Program AANDB Sample Output 
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OATE* 12/01/75 


LIST 

PARAMETER 

table t 

1 1 

1 1 1 

l 1 1 1 1 1 

CASE NUMBER 1 PLOT COOE 

0 

TIME RESPONSE COOE 1 


INPUT 

DATA 




PN, 

TIME. 

0. 000 
0.00 

.820 

.50 

1.020 

1.00 

1.0 20 
1.50 

1*005 1.030 

2. Of 2.50 

BET AN, 
TIME. 

0.000 

0.00 

. 160 
.50 

.580 

1.00 

.9 30 
l. 50 

,980 .800 

2.00 2.50 

OSTAR. 

TIME, 

0. 000 
0.00 

.376 

.50 

.676 

1.00 

.951 

1.50 

1,260 1.610 

2. 0 C 2.50 


1 1 

111 

1 1 



1.070 
3, 00 

1.090 

3.50 

1.075 

4.00 

1.045 

4.50 

1. 025 
5.00 

• 650 
3.00 

.690 

3,50 

. 890 
4.00 

1.090 

4.50 

1.150 
5. no 

1.90 0 
3.00 

2.340 

3.50 

2. 650 
4. 00 

2.950 

4.50 

3.270 

5.00 
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ORIGINAL PAGE IS 
OF POOR QUALITY. 
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REAL 


IMAG 


REAL 


I HAG 


HU, ,15943«E*01 0. 

EX, .9330 0.0000 


• 279846E + 00 0. 

-2.5470 0.0000 


REAL 

. 31766 2E*00 
-,4112 


IHAG 

• 749635E* OO 
2. 3400 


REAL 

,317662 E*00 
-.4112 


IHAG 

• .7 49 635E + 00 
-2.3400 



REAL 

IMAG 

PEAL 

IMA C- 

HU, 

,932945E*00 

0. 

-• 5152 20E^ 01 

0. 

EX, 

-.1388 

0.0000 

3.2788 

6.2832 


REAL 

I MAS 

REAL 

IHAG 

497 955E*00 

• 78 1 25 6E* 0 0 

• 497955E + 00 

-. 781256E*-00 

-.1528 

2.0067 

-.1528 

-2.0067 


real 


IMAG 


PEAL 


IMAG 


MU* . 10 3636E *01 0. ,392694E*00 

El* .0714 0.0000 *1.8694 


REAL 

IMAG 

REAL 

?5 6 11 BE* 00 

• 1O5674E + 0 1 

.25611BEf00 

.1675 

2. 6660 

.1675 
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0.0000 

IMAG 

105674Em 

2.6660 
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0.0 000 
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IMAG 
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• 4515 62E»0Q 

756023£»00 

• 4515 62E* 00 

.756023E4-00 

-.2543 

-2.0648 

-.2543 

2. 0648 


THE GENERATED PHI DATA IS, 
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0. 000 
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TIME, 
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2.50 


2.772 
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5. 00 
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fitting results 
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the a matrix is 


235 877E* 01 
•» 53<#355 E-8i 
«293 87<*e-0t 
• 959S79 E*00 


. 7708322*00 
358690 E *00 
-.100515E*Q1 
-.207575E-01 


-,109 39 OE* C 2 
. 384788E*01 
201315E*tO 
. 276157E* CO 


255620E- 02 
-.2060 77E-03 
.5327 88E-C1 
. 262703E- 0? 


THE 0 VECTOR IS 

. setsoaEm 

-.5273 09E-01 
-,2i5i 79E*00 
.t»95507E-Oi 
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RESPONSE ENVELOPES 


TIME 

0 . 00 

• to 

• 20 

• 30 

• 40 
. SO 

• 60 

• 70 

• 60 
• 90 

1 . 00 

1. tO 

1. 20 
la 3 0 
1*40 

1.50 
1.60 
ie 70 
1*6 0 

1.90 
2*00 

2 . 10 

2. 20 

2. 30 

2.40 

2 . 50 

2. 60 

2.70 

2. BO 

2.90 
3*00 

3. 10 
3.20 

3. 30 

3. 40 

3.50 
3.60 

3.70 

3.60 

3, 90 

4. 00 
4, 10 

4. 20 

4. 30 

4.40 

4.50 

4.60 


PN 6 ETAN DSTAR 


0.000 

0.000 

2.000 

• 333 

• 010 

,000 

.667 

.027 

000 

1.000 

• 065 

2*000 

1.043 

.107 

2.000 

1.064 

• 160 

2.000 

1.076 

• 227 

2,000 

1.060 

• 300 

2.000 

1.060 

• 360 

2.000 

1.078 - 

.409 

2.000 

1.076 

.450 

2.000 

1,074 

.479 

2.000 

1*072 

.507 

2,000 

1.070 

.530 

2.000 

1.069 

.556 

2,000 

1.067 

,579 

2.000 

1.065 

.600 

2.000 

1,063 

.620 

2,000 

1.061 

.639 

2.000 

1.099 

.656 

2.000 

1.057 

.677 

2 . 0 C 0 

1.055 

• 692 

2.100 

1.053 

.70 8 

2,200 

1.051 

.722 

2.300 

1.050 

.733 

2.400 

1.046 

.748 

2.500 

1.046 

. 762 

2.600 

1.044 

. 778 

2 . 700 

1.042 

.790 

2.800 

1 . 040 

. 804 

2.900 

1.038 

.816 

3.000 

1.036 

• B 29 

3.100 

1.034 

.840 

3.200 

1.0 32 

.852 

3.300 

1.030 

. 861 

3.400 

1.029 

.871 

3.500 

1.027 

. 880 

3.600 

1,025 

.688 

3.700 

1.023 

.896 

3.800 

1.021 

. 90 8 

3.900 

1,019 

.915 

4.000 

1.017 

.920 

4.100 

1.015 

. 929 

4 . 200 

1.013 

.93 8 

4.300 

1.011 

.947 

4 o 40 0 

1.010 

.956 

4.500 

1 , 008 

. 965 

4.600 


- 2.0 00 

4.000 

- 4 . 008 

* 2.000 

4 . COO 

- 4.0 00 

- 2.0 00 

4.000 

- 4,000 

- 2 , 000 

4.000 

- 4.0 on 

- 2.0 00 

4.000 

- 4 . 000 

- 2.000 

4.000 

- 4.000 

- 2 . 0 00 

4.000 

- 4.0 00 

- 2. 000 

4 . 000 

- 4 . 000 

- 2.0 00 

4,000 

- 4,000 

- 2 . 0 00 

4 . 000 

- 4. 0 00 

- 2.0 00 

4.000 

- 4. 0 00 

- 2.0 00 

4.000 

- 4.000 

- 2 . 0 00 

4.000 

- 4.000 

- 2 . 0 00 

4.000 

• 4.000 

- 2.000 

4.000 

- 4.0 00 

- 2.000 

4,000 

- 4.000 

- 2 . 000 

4 . 00 C 

- 4 . 0 00 

- 2.000 

4,000 

- 4.000 

- 2 . 0 00 

4 , 000 

- 4 , 0 00 

- 2.0 00 

4,000 

- 4.000 

- 2.000 

4.000 

- 4. 0 00 

- 2.100 

4,200 

- 4.200 

- 2.200 

4,400 

- 4.400 

• 2 . 3 CO 

4.600 

- 4.600 

- 2.4 00 

4,800 

- 4.810 

- 2.500 

5 , 000 

- 5.0 00 

- 2.6 00 

5.200 

- 5,2 00 

- 2.700 

5,400 

- 5.400 

- 2.800 

5,600 

- 5 . 600 

- 2.9 00 

5.800 

- 5.800 

- 3.000 

6,000 

- 6.000 

- 3 . 100 

6,200 

- 6.200 

- 3.200 

6 , 4 0 0 

- 6.400 

- 3.300 

6,600 

- 6.600 

- 3.400 

6,800 

- 6.800 

- 3.5 00 

7.000 

- 7 . 0 00 

- 3 . 600 

7 .20(1 

- 7.200 

- 3 , 700 

7,400 

- 7.400 

- 3.800 

7.600 

- 7 . 6 00 

- 3.9 00 

7.800 

- 7.800 

- 4 . 000 

8,000 

- 8.000 

- 4.100 

8.200 

- 6.200 

- 4 . 200 

8 . 400 

- 8.400 

- 4.3 00 

8.600 

- 8.600 

- 4 . 400 

8 . 800 

- 8.8 00 

- 4.500 

9.000 

- 9 . 000 

- 4,600 

9.200 

- 9.200 
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4 .70 
4.80 
4.90 
5.00 


1 * 006 

• 974 

4 . 700 

- 4 . 7 00 

9,400 

- 9,400 

1*804 

.983 

4.800 

- 4.800 

9.600 

- 9 . 600 

1. 002 

.992 

4.900 

- 4.9 00 

9,800 

- 9.8 00 

t. 000 

1.000 

5,000 

-5, o no 

10,000 

- 10 . 000 


.116 


PND OT 


OSTARDOT 


QtTANOOT 


4.000 

• 006 

6.000 

4.000 

• 711 

6.000 

i . €06 

• 255 

4.800 

t.m 

••224 

4.150 

• 992 

-.297 

3.720 

• 908 

-.323 

3.280 

.643 

-.327 

2.910 

• 797 

-.323 

2.610 

• 224 

-.319 

2.370 

• 663 

-.304 

2 . 100 

• 613 

-.289 

1.950 

.556 

-.274 

1.780 

• 510 

-.251 

1.610 

• 464 

-.228 

1.490 

• 421 

-.205 

1.350 

• 379 

-.183 

1 . 260 

• 349 

-.148 

1.170 

• 318 

-.129 

1.120 

• 295 

-.125 

1.070 

• 272 

-.118 

1.950 

.257 

-.125 

1.800 

• 234 

-.123 

1.000 

• 222 

-.121 

1 . 000 

• 211 

-.120 

1.000 

• 199 

-.ita 

1.000 

• 192 

-.116 

1.000 

• 198 

-•114 

t .000 

• 184 

-.112 

1.000 

.160 

-.110 

1.000 

.176 

-.109 

1.000 

.172 

-.107 

l.OOfl 

• 166 

-.105 

1.000 

.165 

-.103 

1.000 

• 169 

-.102 

1.000 

.165 

-.100 

1.000 

• 165 

-.0 98 

1.000 

.157 

-.096 

1.000 

• 142 

-.095 

1.000 

.119 

-.093 

1.000 

• 107 

-.091 

1.000 

• 096 

-.089 

1.000 

• 060 

-.086 

1.000 

• 073 

-.066 

1.000 

• 069 

-.084 

1.000 

« 061 

-.082 

1.000 

• 054 

-.081 

1.000 

• 046 

-.079 

i.ooo 

• 046 

-.077 

1.000 

. 046 

-.075 

1.000 


- 6.000 

12.000 

- 12.000 

- 6.000 

12.000 

- 12.000 

- 4.800 

9.660 

• 9.600 

- 4.150 

8.300 

- 8,3 00 

- 3.720 

7.440 

- 7.440 

- 3.280 

6.560 

- 6.5 60 

- 2.910 

5 * 820 

- 5.820 

- 2.610 

5.220 

- 5.220 

- 2.370 

4,740 

- 4.7 40 

- 2.100 

4*200 

- 4.2 00 

- 1.950 

3.900 

- 3.9 00 

- 1.780 

3.560 

- 5.560 

- 1 s 610 

3*220 

- 3.220 

- 1.490 

2.960 

- 2.980 

- 1.350 

2.700 

- 2.7 00 

- 1.260 

2.520 

- 2.5 20 

- 1.170 

2*340 

- 2.340 

- 1.120 

2.240 

- 2.240 

- 1.070 

2.140 

• 2.140 

- 1.050 

2.100 

• 2.100 

- 1.000 

2.000 

- 2.0 00 

- 1.000 

2.000 

- 2,000 

- 1.000 

2.000 

- 2.0 00 

- 1.000 

2.000 

- 2 , 00(1 

- 1.000 

2.000 

- 2.0 00 

- 1.000 

2.000 

- 2.0 00 

- 1.000 

2.000 

- 2.000 

• 1.000 

2.000 

- 2.000 

- 1.000 

2. 0 0 i 3 

- 2.000 

- 1.000 

2.000 

- 2.000 

- 1.000 

2.000 

- 2.0 00 

- 1.000 

2.000 

- 2.000 

- 1.000 

2.000 

- 2.000 

-1 .000 

2.000 

- 2.000 

- 1.000 

2,0 00 

- 2.000 

- 1.000 

2.000 

- 2.000 

- 1.000 

2.000 

- 2.000 

- 1.000 

2.000 

- 2.000 

- 1.000 

2.000 

- 2,0 00 

- 1.000 

2.000 

- 2.000 

• 1.000 

2.000 

- 2.000 

- 1.000 

2.000 

- 2.000 

- 1.000 

2.000 

- 2,0 00 

-1 .000 

2.000 

- 2,000 

- 1.000 

2.000 

- 2,0 00 

- 1.000 

2.000 

- 2,000 

- 1.000 

2.000 

- 2.000 

- 1.000 

2.000 

- 2.000 

- 1.000 

2.000 

- 2.000 
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* 046 

• Sir6 


••074 1.000 - 1.000 2.000 

*.072 1.000 - 1.000 2.000 


- 2.000 
- 2.0 00 
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fitteo time responses 


TIME 

PN 

0 E T AN 

0 . 00 

• 000 

-.000 

• 10 

• 252 

. 003 

.20 

• 450 

r 020 

• 30 

• 606 

.0 52 

• 40 

.727 

• 099 

• 50 

• 819 

. 159 

• 60 

.889 

.231 

.70 

• 941 

.312 

• 80 

• 978 

.399 

• 90 

1.003 

. 4 J 9 

1.00 

1 • 01 B 

.5 79 

1.10 

1.027 

.666 

i . 20 

1.030 

. 747 

1.30 

1.029 

.819 

1 . 40 

1*026 

. 881 

1.50 

1*021 

.9 30 

1.60 

1.017 

.966 

i . 70 

1.012 

. 989 

1.80 

1.009 

.998 

1 « 90 

1.006 

. 995 

2.00 

1.006 

.980 

2 . 10 

1.007 

. 955 

2.20 

1.010 

. 923 

2 . 30 

1.015 

.885 

2.40 

1.021 

. 844 

2.50 

1.028 

.802 

2.60 

1.037 

.761 

2.70 

1.045 

.723 

2 . 80 

1.054 

.6 91 

2.90 

1.063 

. 665 

3 . 00 

1.070 

.648 

3. 10 

1.077 

.63 9 

3 . 20 

1.083 

.639 

3 . 30 

1.087 

.648 

3.40 

1.090 

. 666 

3 , 50 

1.092 

.691 

3 . 6 G 

1.091 

.724 

3.70 

1.090 

.762 

3.80 

1 . 086 

. 80 4 

3.90 

1.082 

• 8 48 

4 . 00 

1.077 

. 894 

4 . 10 

1.071 

. 939 

4 . 20 

1.064 

. 982 

4.3 0 

1.057 

1 . 022 

4.40 

1. 050 

1.057 

4.5 0 

1 . 044 

1.0 88 

4.60 

1.038 

1 . 112 


OSTAR 

PNOOT 

9ETAN00T 

-.000 

2.828 

-.038 

. 082 

2.232 

.0 97 

.160 

1.754 

.245 

.236 

1.369 

.395 

. 308 

1. 058 

.538 

. 376 

• 80 6 

.665 

• 441 

.601 

. 769 

.504 

.436 

.846 

. 563 

.304 

. 893 

.621 

.199 

.907 

.677 

.117 

. 690 

.732 

.055 

• 843 

.786 

. 009 

. 769 

• 640 

-.021 

. 671 

.895 

-.040 

.556 

. 951 

-.048 

.428 

1.0 09 

-.048 

. 294 

1 . P68 

-.0 41 

• 158 

1 . 129 

-.029 

.028 

1.192 

-.013 

-.09 3 

1.257 

. 004 

-. 199 

1.325 

• Q22 

-.288 

1 . 394 

.039 

-.355 

1.465 

.055 

-. 400 

1.538 

.068 

-. 421 

1.611 

.070 

419 

1.686 

.085 

-. 396 

1.761 

• 088 

-. 352 

1 .835 

. 087 

-. 29 2 

1. 910 

.0 82 

-.218 

1.983 

.075 

-. 134 

2.056 

.064 

-. 045 

2 . 128 

.051 

. 046 

2.198 

.037 

. 135 

2.267 

. 021 

. 2iB 

2.334 

.005 

.292 

2 . 400 

-.011 

. 355 

2.464 

-.0 25 

. 403 

2.527 

-.0 38 

. 437 

2.590 

-.0 49 

. 454 

2.651 

-.050 

. 456 

2.712 

-.064 

. 443 

2.772 

-.067 

• 416 

2 . 832 

-.068 

. 377 

2.893 

-.066 

. 329 

2.953 

-.062 

.273 

3.015 

-.055 

. 21 3 


OSTAROOT 

.832 

♦ 803 
. 7*0 
.736 

♦ 70 2 

♦ 66 6 

♦ 637 
.609 
.565 
. 566 
.553 

• 545 
.542 
. 545 
.553 
o 566 
. 532 
.601 
.621 
.643 
.664 
.68 4 
.703 
. 719 
.731 
. 741 
.746 
.748 
.746 
.741 
. 733 
.722 
.709 
.695 
.630 
.665 
.651 
.63 8 
.626 
.617 
.609 
.605 
.603 
.603 
. 60 6 
.611 
.617 
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WO 1*033 

<*•00 1,029 

<*.90 1.026 

0.00 1.024 


1 . 130 3.077 

1.142 3.140 

1.148 3.204 

1 . 149 3.269 


-.047 • 151 

-.037 .090 

-.027 . 032 

-.016 -.019 


.635 

.634 

.644 

.634 
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INTEGRATED RESPONSE 


TIME 

PN 

o* oo 

0.000 

• 10 

.252 

• 20 

.450 

• 30 

*6 06 

• 40 

.727 

*50 

• 8 19 

*60 

• 889 

*ro 

.9 41 

• 80 

.978 

• 90 

1.003 

1*00 

1.018 

1 . 10 

1.027 

1.20 

1 . C 30 

1 . 30 

1.0 29 

1.40 

1.0 26 

1.53 

1.0 21 

1.60 

1.017 

1*70 

1.012 

1.80 

1.0 09 

1.90 

1.0 06 

2 * 00 

1.0 06 

2 . IQ 

1.0 07 

2*20 

1.010 

2.30 

1.015 

2 . 40 

1.021 

2 . 50 

1.028 

2.60 

1.037 

2*70 

1.0 45 

2.80 

1.054 

2.90 

1.0 63 

3 . 00 

1.070 

3 . 10 

1.077 

3.20 

1.083 

3 . 30 

i.oar 

3.40 

1.090 

3 . 50 

1.0 92 

3.60 

1.091 

3.70 

1.090 

3.80 

1.086 

3.90 

1.082 

4 . 00 

1.0 77 

4 * 10 

1.071 

4.20 

1.0 54 

4.30 

1.057 

4.40 

1.0 50 

4.50 

1 . 044 

4 . 60 

1.038 


SET AN 

DSTAR 

0.000 

-.0 00 

• 003 

• 082 

.020 

• 160 

.0 52 

• 236 

.099 

* 3 08 

.159 

.176 

.231 

.441 

.312 

.504 

.399 

.563 

*489 

.521 

.5 79 

.677 

. 666 

.732 

.747 

• 7 86 

.819 

.840 

.881 

. 895 

,930 

.951 

.966 

1.009 

.989 

1.068 

.998 

1. 129 

.995 

1.192 

.980 

1.257 

.955 

1.325 

.923 

1.3 94 

• 88 5 

1.465 

.844 

1.538 

.802 

.1. 611 

.761 

1.686 

.723 

t. 761 

.691 

1.835 

.665 

1. 910 

.648 

1.983 

.639 

2.056 

.639 

2. 128 

.648 

2.198 

.666 

2.267 

.691 

2.334 

.724 

2.40Q 

.762 

2.464 

.804 

2.527 

.848 

2.590 

.894 

2.651 

.9 39 

2.712 

.982 

2.772 

1.0 22 

2. 8 32 

1.057 

2.893 

1.0 8-8 

2.953 

1.112 

3.015 


4 . 70 
4.60 
4 . 90 
6.00 


1.033 
1.0 29 
1 . 026 
1.0 24 


1.130 

1.142 

1.148 

1.149 


3.0 77 
3.140 
3.204 
3 . 269 
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FIRST DERIVATIVES OF 


TIME 

PNOOT 

0.03 

2,928 

• 10 

2.2 32 

• 20 

1.754 

• 30 

1.369 

• 40 

1.0 58 

• 50 

.9 06 

• 60 

• 601 

• 70 

• 436 

• eo 

.3 04 

• 90 

• 199 

1. 00 

.117 

1* 10 

.0 55 

1.20 

.0 09 

1, 3(1 

-.0 21 

1.40 

-.0 40 

1, 50 

-.0 40 

1.60 

-.048 

l.ro 

-*041 

i.ao 

-.029 

1.90 

-.0 13 

2 , 00 

• 004 

2.10 

• 0 22 

2.20 

.039 

?• 30 

.055 

2*40 

.0 68 

2.50 

.078 

2.60 

.085 

2.70 

.0 86 

2. BO 

. 087 

2. 90 

• 0 32 

3 a «0 

.0 75 

3.10 

.0 64 

3. 20 

• 0 51 

3. 30 

.037 

3,40 

.021 
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